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ABSTRACT 


It  has  been  shown  that  good  correlation  exists  between  large-scale,  high  level 
signal  returns  of  bottom  acoustic  reverberation  and  bathymetric  ridge  structures.  The 
ultimate  goal  of  this  type  of  analysis  would  be  the  removal  of  propagation  effects 
resulting  in  a  large-scale  mapping  of  scattering  strengths.  Furthermore,  analysis  with  a 
quasi-CW  propagation  model  suggests  a  strong  correlation  between  the  small-scale 
fluctuations  in  the  reverberation  signal  and  the  bottom  acoustic  ensoniflcation,  thereby 
suggesting  that  uncertainties  in  the  predictions  of  forward  propagation  may  limit  the 
resolution  of  such  mapping.  By  employing  broadband  modeling  techniques,  a  valid 
representation  of  the  complete  time  domain  forward  propagation  is  provided.  Diffuse 
secondary  bottom  interactions  appear  to  affect  the  resolution  of  the  primary,  direct-path 
interaction  at  ranges  beyond  a  few  water  depths.  Analysis  of  data  recorded  by  near¬ 
bottom  vertical  line  arrays  (VLA’s)  confirms  the  existence  of  these  secondary,  multipath 
interactions  in  the  forward  propagation.  The  exact,  two-way  travel  times  fi-om  all  bottom 
interactions  are  modeled  and  the  influence  of  the  secondary  interactions  is  quantified. 
Possible  ramifications  for  general  sonar  system  performance  are  discussed. 
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I.  INTRODUCTION 


In  recent  research  conducted  using  data  collected  under  the  Acoustic  Reverberation 
Special  Research  Project  (ARSRP)  sponsored  by  the  Office  of  Naval  Research,  it  has 
been  shown  that  good  correlation  exists  between  large-scale,  high  level  signal  returns  of 
bottom  acoustic  reverberation  and  bathymetric  ridge  structures.  One  may  then  attempt 
to  obtain  measures  of  bottom  scattering  strengths  by  removing  the  influence  of 
propagation.  Analysis  with  a  quasi-CW  propagation  model,  however,  suggests  a  strong 
correlation  between  the  small-scale  fluctuations  in  the  reverberation  signal  and  the  bottom 
acoustic  ensonification.  This  suggests  that  uncertainties  in  the  predictions  of  forward 
propagation  may  limit  the  resolution  of  such  an  inversion.  The  limits  of  the  resolution 
of  such  a  map  are  also  imposed,  in  part,  by  the  processing  and  the  knowledge  of  the 
environment. 

By  employing  broadband  modeling  techniques,  a  valid  representation  of  the 
complete  time  domain  forward  propagation  can  be  obtained.  Despite  the  use  of  a 
directional  vertical  source  array,  diffuse  secondary  bottom  interactions  may  affect  the 
resolution  of  the  primary,  direct-path  interaction  at  ranges  beyond  a  few  water  depths. 
Analysis  of  data  recorded  by  near-bottom  vertical  line  arrays  (VLA’s)  confirms  the 
existence  of  these  secondary,  multipath  interactions  in  the  forward  propagation.  By 
computing  the  complex  field  at  the  bottom  interface  from  both  a  vertical  line  array  of 
sources  and  a  point  receiver,  co-located  in  a  monostatic  geometry,  the  exact,  two-way 
travel  times  and  transmission  loss  from  all  bottom  interactions  can  be  determined.  In  this 
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manner,  the  influence  of  the  secondary  interactions  on  the  measured  bottom  reverberation 
can  be  estimated. 

The  purpose  of  this  thesis  is  to  examine  these  effects  and  attempt  to  quantify  them. 

Other  issues  are  what  level  of  resolution  exists  within  the  measured  reverberation  at 
ranges  around  one-half  of  the  convergence  zone  distance  and  greater,  and  what  level  of 
resolution  can  model  predictions  hope  to  provide.  Schemes  for  removing  the  influence 
of  these  secondary  multipath  interactions,  if  necessary,  will  also  be  discussed.  Addressing 
these  issues  will  advance  the  general  knowledge  of  underwater  acoustic  propagation. 
Quantification  of  these  influences  and  suggesting  schemes  to  overcome  them  will  lead  to 
future  improvements  in  sonar  system  design.  The  ability  to  resolve  complex  multipaths 
at  significant  ranges  in  a  covert  maimer  would  provide  great  advantages  to  units  involved 
in  surveillance/power  projection  operations. 

In  Chapter  H,  some  background  information  on  the  ARSRP  experiment  is 
provided.  The  numerical  model  used  to  predict  acoustic  propagation  is  also  described 
along  with  other  processing  techniques  used  in  the  data  analysis.  In  Chapter  III,  the 
existence  of  the  secondary,  multiple  arrivals  is  confirmed  in  the  measured  data  and 
compared  with  the  forward  propagation  model  predictions.  The  exact,  two-way  travel 
times  and  transmission  loss  are  then  modeled  and  the  influence  of  the  secondary 
interactions  is  quantified.  Findings  are  summarized  in  Chapter  IV  and  possible 
ramifications  for  general  sonar  system  performance  are  discussed. 
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11.  BACKGROUND  INFORMATION 


A.  THE  ARSRP  STUDY 

In  1990,  the  Office  of  Naval  Research  (ONR)  began  an  Acoustic  Reverberation 
Special  Research  Program  (ARSRP)  in  the  area  of  underwater  acoustic  reverberation.  The 
scientific  plan  for  the  bottom  reverberation  component  of  the  ARSRP  placed  special 
emphasis  on  obtaining  a  detailed  description  of  the  bottom  in  the  "natural  laboratory" 
where  acoustic  experiments  were  to  be  conducted  [Ref.  1].  The  location  of  the  ARSRP 
natural  laboratory  was  selected  to  be  on  the  west  flank  of  the  Mid- Atlantic  Ridge 
(approximately  25"  to  27"  N  and  45"  to  50"  W).  Adequate  environmental  measurements 
could  be  obtained  here  in  order  to  understand  the  connection  between  sea  floor 
characteristics  and  scattering  characteristics.  Two  geology  and  geophysics  surveys  were 
conducted  of  the  region  which  resulted  in  data  gridded  at  200m  ^  200m.  More  detailed 
analysis  was  conducted  at  locations  of  particular  interest  providing  information  on  the 
small-scale  roughness  spectra  and  bottom  composition. 

The  ARSRP  has  conducted  two  acoustic  field  experiments  in  the  region.  An 
acoustic  reconnaissance  experiment  was  carried  out  in  July- August  1991  aboard  the 
Research  Vessel  (R/V)  Cory  Chouest  [Ref.  2].  The  major  bottom  reverberation 
experiment  of  the  ARSRP  was  carried  out  in  July  1993  [Ref.  3].  Three  ships  participated 
in  the  experiment  -  the  R/V  Cory  Chouest,  the  R/V  Alliance,  and  the  R/V  Knorr.  The 
Cory  deployed  a  vertical  line  source  array  and  a  horizontal  line  receive  array,  the 
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Alliance  deployed  a  flextensional  source  (vertical  pair)  and  a  horizontal  line  receive  array, 
and  the  Knorr  deployed  both  near-bottom  vertical  line  receive  arrays  and  recording 
packages  as  well  as  the;  DTAGS  (Deeply  Towed  Acoustics  Geophysics  System).  Fig.  1 
shows  the  ship  tracks  for  the  Cory  and  Alliance.  A  major  emphasis  of  the  experiment 
was  on  the  collection  of  bistatic  scatter  from  Sites  A,  B’,  and  C’.  Both  the  ship  tracks  and 
the  transmission  sequences  for  the  Cory  and  Alliance  were  designed  so  that  each  ship 
would  collect  monostatic  scatter  from  their  own  pings  along  with  bistatic  scatter  from 
each  other. 

During  the  1993  acoustics  experiment,  both  the  Cory  and  the  Alhance  made  water 
colunrn  temperature  and  conductivity  measurements.  The  Cory  deployed  63  expendable 
bathythermographs  (XBTs)  as  well  as  four  expendable  CTDs  (1000m).  The  Alliance 
deployed  39  XBTs  and  made  conventional  CTD  casts  to  full  ocean  depth  at  five  locations. 

B.  THE  BEAMFORMER 

TTie  beamformer  developed  to  process  the  transmission  loss  data  output  from  the 
acoustic  propagation  model  was  written  as  a  Matlab  m-file.  The  m-file  has  been  written 
to  process  the  data  from  a  63-element  vertical  line  array,  which  simulates  the  near-bottom 
vertical  arrays  used  in  the  ARSRP  study.  Inter-element  spacing  of  the  simulated  array 
is  1.95  meters.  The  actual  arrays  used  in  the  study  consisted  of  59  elements,  with  an 
inter-element  spacing  of  2.07  meters.  The  important  consideration  for  comparison  of 
simulated  and  actual  data  is  that  their  positions  and  total  lengths  are  in  agreement.  The 
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Figure  1.  Ship  tracks  of  R/V’s  Cory  Chouest  and  Alliance  during  1993  main  acoustics  experiment.  From  Ref.  [4]. 


data  is  processed  to  produce  a  vertical  arrival-angle  versus  time  magnitude  plot  of  the 
received  sonar  signal  at  the  vertical  array. 

The  beamformer  is  written  for  the  specific  output  of  the  propagation  model’s 
broadband  simulations.  This  data  corresponds  to  a  frequency  bandwidth  of  66.66  Hz  with 
a  center  frequency  of  227.5  Hz  and  a  total  of  1024  frequencies.  Thus,  incremental 
frequency  spacing  is  66.66/1024  =  0.065104  Hz.  This  results  in  a  total  Fourier- 
transformed  time  window  of  15.36  seconds  and  a  time  resolution  of  0.015  seconds. 

These  1024  frequencies  are  placed  in  the  positive  portion  of  a  frequency  window 
which  extends  from  -533.33  Hz  to  +533.33  Hz,  with  all  other  frequencies  zero  padded. 
A  Blackman  window  is  applied  to  the  non-zero  frequency  information,  and  then  the  entire 
frequency  window  is  transformed  to  the  time  domain  resulting  in  a  time  window  of  15.36 
seconds  long,  but  a  time  resolution  of  9.4  x  10"^  seconds. 

This  time  window  is  broken  into  thirty-one  overlapping  windows  to  take  advantage 
of  the  increased  time  resolution  and  minimize  sidelobe  effects  m  time.  Each  overlapping 
window  consists  of  1024  sampled  times,  which  are  transformed  back  to  the  frequency 
domain  where  the  frequency  resolution  is  now  1.043  Hz.  A  Blackman  window  is  now 
applied  to  the  63-element  spatial  array  and  transformed  to  the  vertical  wavenumber 
domain.  Arrival  angle  information  corresponds  to  wavenumber  information  as  follows 
[Ref.  5]: 


m  = 


(N+Z)dcost|j 

m 

I 


(N+Z)dcos\|r  f 

m 

C 


(1) 
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where  m  is  the  wavenumber  bin,  N  is  the  number  of  array  elements,  Z  is  the  number  of 
zeros  that  the  FFT  is  padded  with,  d  is  the  inter-element  spacing,  i|rn,  is  the  arrival  angle 
relative  to  endfire,  A  is  the  wavelength,  f  is  the  frequency,  and  c^,  is  the  reference  sound 
speed  of  the  water. 

Once  the  arrival  angle  information  is  determined,  the  frequency  window  is  finally 
transformed  back  to  the  time  domain.  The  time  domain  information  in  the  overlapping 
portions  of  the  time  windows  are  incoherently  summed  and  the  final  matrix  of  arrival 
angle  versus  time  is  converted  to  transmission  loss  in  dB  (unnormalized)  for  plotting. 

C.  THE  ACOUSTIC  PROPAGATION  MODEL 


The  full-wave  propagation  method  used  to  model  the  acoustic  propagation  for 
these  studies  uses  a  parabolic  equation.  Since  the  introduction  of  the  parabolic  equation 
method  into  the  field  of  underwater  acoustics  in  the  1970’s  by  Hardin  and  Tappert  [Ref. 
6],  the  PE  method  has  become  a  popular  wave-theory  technique  for  solving  range- 
dependent  propagation  problems  in  ocean  acoustics. 

1.  Standard  PE  Derivation 

The  starting  point  in  the  derivation  of  the  standard  parabolic  equation  is  the  free- 
field  Helmholtz  equation  for  a  constant-density  medium  in  cylindrical  coordinates  (r,(|),z) 
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and  for  a  harmonic  field  of  time  dependence  exp(-io)t). 


a^p  lap  a^p  2  2  „ 

—  + - + - +k  n  p=  0, 

2  r  ar  2  0  ^  ’ 


ar 


dz‘ 


(2) 


where  azimuthal  symmetry  is  assumed,  p(r,z)  is  the  acoustic  pressure,  =  co/Cq  is  a 
reference  wavenumber,  n(r,z)  =  C(/c(r,z)  is  the  index  of  refraction,  and  c(r,z)  is  the  sound 
speed  describing  the  environment. 

Assume  the  solution  to  Eq.(2)  takes  the  form 


p(r,z)=  r) , 


(3) 


which  is  an  outgoing  cylindrical  wave  solution.  The  envelope  function  i}r(r,z)  is  assumed 
to  be  slowly  varying  in  range.  The  Hankel  function,  which  satisfies  the  Bessel 
differential  equation 


a^iy^”(rr)  .  dH^'\k  r) 

0  0^  1  0^0'^ 


(1). 


dr^ 


dr 


^klH^\kx)=  0, 

0  0  0 


(4) 


can  be  replaced  by  its  asymptotic  form  for  )fe(,r  »  1, 


%  r)  = 


Tzk  r 
0 


i(*^r-ir/4) 

e 


(5) 
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Substituting  Eq.(3)  into  Eq.(2),  and  using  the  Hankel-function  property  of  Eq.(4),  gives 


H^\k  r) 
0  ^  0 


ar 


1^ 
+  — 

ry 


aiif  a^ijf  ,  2  2 

—  +  — -+k  (n  -1)1];  =  0 .  (6) 

^  dz^  ® 


Making  the  far-field  assumption,  »  1,  Eq.(5)  can  be  used  to  obtain  the  simplified 
elliptic  wave  equation 


a^ijt 

ar^ 


atjt 

+  Tik  —  + 

0  ar 


a^i|t 

dz^ 


(7) 


Next,  the  paraxial  approximation  is  used  to  arrive  at  the  standard  parabolic  wave 
equation.  This  small-angle  approximation  is  expressed  by 


a^t|j 

ar^ 


aijr 

clik  — . 

0  ar 


(8) 


A  qualitative  justification  for  this  approximation  can  be  made  by  noting  that  the  main 
radial  dependence  of  the  field  is  contained  in  the  Hankel  function  through  the  term 
exp(i^),  while  the  envelope  tj;  will  vary  slowly  with  range  over  a  wavelength  A.  This 
can  be  expressed  mathematically  as  0t|f/ar  <  ~  ijr ,  or  a^x|f/ar^<  ik^dT}f/dr ,  and 

consequently  the  paraxial  approximation  is  justified. 
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Making  use  of  Eq.(8)  in  Eq.(7)  results  in 


2ik 


0ilr 

di 


dz^ 


*  (n"-l)l|l  =  0 , 


or 


(9) 


3t|r 

"aT 


i  a^tjj 


ik 


which  is  the  standard  parabolic  equation. 


2.  The  Split-Step  Fourier  Algorithm 


A  second-order  differential  equation  in  r  has  now  been  reduced  to  a  first-order  one 
in  r,  allowing  solutions  by  way  of  a  non-iterative  marching  algorithm.  Eq.(9)  may  be  re¬ 
written  as 


i  aijr 

k  di 
0 


t  =  (r 


op 


op 


(10) 


where  the  operators 


T 

op 


i  d 


\2 


k  dzj 
0 


1 

2k^  dz^ 
0 


(11) 
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and 


^op  "  =  --i  (n^  -1) .  (12) 

This  representation  of  the  operators  as  kinetic  and  potential  energy  operators  is 
especially  insightful  when  it  is  desired  to  form  the  ray  equations  which  have  Hamiltonian 
form  [Ref.  7].  In  Eq.(lO),  the  function  ij;  is  a  vector  (in  z)  in  Hilbert  space.  The  values 
of  t|;  can  now  be  expressed  as 

i|f(r  +  Ar)  =  $(r)tlr(r).  (13) 

For  all  propagation  predictions  presented,  we  have  employed  the  University  of 
Miami  Parabolic  Equation  (UMPE)  model  of  Smith  and  Tappert  [Ref.  8].  The  UMPE 
model  uses  the  spUt-step  Fourier  (SSF)  method  [Ref.  6]  to  compute  the  PE  solution.  This 
is  accomplished  by  approximating  the  propagator  function  ^(r)  by 


where 


$(r)  »  e 


-ik  H  (r)Ar 

0  op^  ' 


H 

op 


1 

Ar 


I 


r  +  Ar 


op 


(14) 


(15) 


This  can  be  approximated  by 


H 

op 


(r  +  i  Ar) , 


(16) 
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or  simply 


(17) 

which  are  referred  to  as  "centered"  and  "endpoint"  schemes,  respectively.  Presumably  if 
the  range  step  Ar  is  small  enough  the  differences  between  the  solutions  are  negligible. 
The  UMPE  model  uses  the  approximation  of  Eq.(17). 

The  operator  U^p  is  a  multiplication  operator  in  z-space  and  is  therefore  a  diagonal 
matrix.  The  operator  T^p  is  not  diagonal  in  z-space,  but  is  diagonal  in  wavenumber  space. 
By  separation  of  H^p  into  the  two  components  T^p  and  U^p  the  propagator  function 
becomes 


0(r)  =  e 


-ikAtT  -ik  AiU 
0  ^  ^ 


where  the  Baker-Campbell-Hausdorff  expansion  [Ref.  9] 


(18) 


e  ^  =  e  ^  ^  ^  +  •  •  •  (19) 

has  been  used.  Since  both  T^p  and  U^p  are  small,  their  products  are  of  second  order  and 
assumed  negligible. 

The  general  algorithm  behind  the  PE/SSF  implementation  is  then  as  follows.  The 
PE  field  function  i|f  is  specified  at  some  range  r  in  the  z-domain.  A  transformation  is 
made  to  the  A:-domain  followed  by  a  multiplication  of  the  k-space  operator  e  . 

The  result  is  then  transformed  again  to  the  z-domain  and  is  followed  by  a  multiplication 
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-ik  LtU 

of  the  z-space  operator  e  .  The  final  result  is  the  field  function  at  r  +  Ar.  The 

FFT  subroutine  employed  in  the  numerical  code  assumes  the  convention 


i|r(z)  =  FFT(tIr(^))  = 


(20) 


and 


m  = 


PFT(i|f  *(z))  *  =  I *  dz . 


Therefore,  the  PE/SSF  implementation  can  be  represented  by 


(21) 


\|f  (r  +  Ar)  =  e 


-ik  AtU  (t,z) 
0 


^  FP 11  e 


X  [FFr(i|(‘(r,z))]‘ 


(22) 


where,  in  fc-space,  the  operator  in  Eq.  (11)  is  given  by 


t  (k)=ll^]\ 

op  l^k' 


(23) 


and  Uop  is  as  defined  in  Eq.  (12). 

The  UMPE  model  allows  the  user  to  choose  from  five  different  sets  of  operator 
forms.  For  the  modeling  used  in  this  thesis,  the  higher  order  wide-angle  approximation 
(WAPE)  form  [Ref.  9]  has  been  used  which  is  accurate  to  greater  than  ±  20°  from  the 
horizontal.  The  operator  forms  for  the  wide-angle  approximation  are 

dz^  ^ 


WAPE 


0 

2 
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and 


U 

WAPE 


=  -(n-1). 


In  wavenumber  space  Eq.(24)  may  be  written  as 


(25) 


-1 


=  1- 


Note  that  modes  with  k  >  are  evanescent  since 


(26) 


(27) 


In  the  modeling  done  for  this  thesis,  measured  sound  speed  profiles  from  the 
Alliance  CTD  casts  are  used  as  input.  The  bottom  is  typically  modeled  as  a  single 
interface  at  depths  extracted  from  the  gridded  database.  The  acoustic  bottom  properties 
are  chosen  to  be  representative  of  this  area  and  have  no  lateral  variability  due  to  lack  of 
geophysical  data.  Various  analytical  formulations  of  source  functions  are  available;  the 
point  source  and  vertical  line  array  were  used  in  this  work.  Broadband  UMPE  model 
runs  were  performed  to  allow  for  more  accurate  investigations  of  the  time  domain 
response  of  the  medium. 
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3.  Computing  One-Way  (Forward)  Travel  Times 


The  UMPE  code  calculates  the  time  domain  arrival  structure  at  a  given  range  in 
a  straight  forward  manner.  The  time  harmonic  acoustic  field  assumed  by  the  model  is 
given  by 


P(r,z,(»)t)  =p  (r,z)e 

(0 


(28) 


The  representation  in  the  time  domain  is  then  simply 


P(r,z,t)  =FFTr[p  (r,z)l=  fp  (r,z)e"'“^do), 

,  u)  J  J  <0 


(29) 


which  follows  from  the  realization  that  Eq.(28)  is  a  single  firequency  component  of  the 
full-spectrum  time  dependent  field  given  by  Eq.(29). 

To  compute  the  arrival  structure  at  some  fixed  range  r=R,  the  model  computes 
and  stores  the  complex  field  values  Pujr(z)  for  many  frequencies  and  then  Fourier 
synthesizes  them  to  obtain  P  (z,t),  the  set  of  complex  pressure  values  in  time/depth 
space.  Because  the  EFT  assumes  inputs  over  the  frequency  band  4  -  BW/2  to  f^  +  BW/2, 
where  BW  is  the  bandwidth  of  the  acoustic  source  and  f^  is  the  center  frequency,  high 
frequency  calculations  can  become  computationally  burdensome  even  for  small 
bandwidth  fields.  Therefore,  the  UMPE  model  frequency  shifts  the  center  frequency  to 
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d.c.  The  calculation  becomes 


_  /  r  -j(u-w  )t 

Pjfet)  =  I  p^^(z)  e  •  d(0,  (30) 

where  o)o  =  27tfo. 

In  practice,  the  model  computes  the  FFT  of  (1/y^ )  ^^(z).  Because  this  neglects 

ik  R  ) 

the  overall  phase  factor  g  °  =  e  * ,  the  time  domain  is  heterodyned  around  the  value 
to  =  R/Co-  Arrival  times  are  then  given  as  values  of  "reduced  time",  or  (t  -  O. 

To  minimize  the  effects  of  sidelobes  due  to  a  flat  frequency  spectrum,  and  of 
aliasing  or  wrap-around  due  to  the  finite  BW,  the  UMPE  model  creates  a  Hamming 
window  which  is  applied  to  the  entire  bandwidth.  This  method  of  tapering  greatly 
reduces  the  effect  of  side  lobes.  Wrap-around  effects  may  still  be  evident  if  the  time 
window  is  not  wide  enough.  The  total  length  of  the  time  window  is 


T  = 


N 

BW 


(31) 


where  BW  is  the  bandwidth  of  the  source  and  N  is  the  number  of  frequencies  for  which 
the  model  computes  solutions,  or  the  frequency  FFT  size  if  the  FFT  is  padded  with  zeros. 
The  frequency  resolution  is 


Ao) 


BW  _  1 

~W~T 


(32) 


so  the  wider  the  time  window  needed,  the  finer  the  bandwidth  must  be  resolved  and  the 
more  total  runs  at  different  frequencies  must  be  performed. 
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4.  Computing  Exact,  Two-Way  Travel  Times 

The  time  domain  arrival  structure  of  the  acoustic  field  at  a  given  point  in  range 
and  depth  can  be  constructed  from  the  Fourier  synthesis  of  the  full-spectrum  complex 
field  values  calculated  by  UMPE  at  that  range  and  depth  over  the  frequency  bandwidth. 

In  this  marmer,  for  every  point  on  the  bottom,  a  time  series  can  be  constructed 
which  would  show  the  arrival  structure  of  the  acoustic  field  from  some  source  at  a 
distance  r  =  R  from  the  point  on  the  bottom. 

Imagine  a  second  source  co-located  at  range  R  from  the  point  on  the  bottom.  If 
we  label  the  bottom  point  b,  and  the  two  co-located  sources  r  and  s,  then  the  arrangement 
would  be  as  shown  in  Fig.  2. 


S 


Figure  2.  Source  and  Receiver  Arrangement 
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The  time  series  at  a  bottom  point  b  from  a  source  at  location  s  could  be 
constructed  as  previously  shown.  Also,  the  time  series  at  receiver  location  r  from  a 
source  at  location  b  could  be  constructed  in  the  same  manner.  By  invoking  reciprocity, 
the  complex  field  at  location  r  from  a  point  source  at  b  over  a  given  frequency  spectrum 
is  the  same  as  the  complex  field  at  location  b  from  the  same  point  source  at  r. 

ITie  time  domain  arrival  stracture  of  the  complex  acoustic  field  at  location  r  from 
a  source  at  point  s  would  be  a  convolution  of  the  time-domain  arrival  structures  from  s 
to  b  and  from  b  to  r: 


=  /ilr^^(t-T,R)  tjr  Jt,R) dx  (33) 


or,  in  the  frequency  domain: 


t^Cco.R)  =  iIr^/w,R)  X  tjr  Ja),R) 


(34) 


where  the  subscript  "2"  stands  for  two-way,  and  t{r^  ((«),R)  =  tjr  ^(a),R)by  reciprocity, 
tir  (&),R)and  tk  (<«),R)  are  the  complex  acoustic  fields  output  by  the  UMPE  model  for 
all  frequencies  in  the  spectrum. 

This  development  has  assumed  reflection  of  tbe  field  at  only  one  bottom  point. 
In  reality,  the  arrival  structure  at  point  r  is  the  summation  of  the  convolved  arrival 
structures  for  all  range  steps: 


2Ar  ,  4Ar 

- ,r+Ar)  +  il;2(t+ - ,r+2Ar)  +  ... 

c  c 

0  o 


(35) 
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where  all  time  series  developed  for  each  step  in  range  must  be  summed.  For  any  given 
range  R,  the  center  time  from  the  Fourier  synthesis  will  be  T,.  =  R/Co  with  a  time  window 
width  of  AT  =  l/6f  where  6 f  is  the  incremental  frequency  spacing  between  each  CW 
calculation.  A  range  step  was  used  to  produce  changes  in  equal  to  an  integer  number 
of  6t  =  1/BW,  the  corresponding  incremental  sampling  of  the  time  window,  i.e. 


In  order  to  develop  the  full  time  series,  the  exponential  factor  immediately  above 
will  be  left  out  of  the  actual  calculation,  and  each  subsequent  range  step  contribution  will 
simply  be  shifted  in  time  so  that  its  center  time  is  t  =  2Ar/c. 
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III.  DATA  ANALYSIS  AND  MODEL  PREDICTIONS 


A.  THE  MODELED  STUDIES  CONDUCTED 


In  order  to  perform  model  predictions  which  can  be  compared  to  actual  five 
segments  were  chosen  from  the  ARSRP  study  [Ref.  3]  which  correspond  to  five  distinct 
arrangements  of  towed  source  transmitting  to  a  near-bottom  vertical  array.  For  all  cases, 
the  model  source  was  placed  at  181  meters  depth,  a  source  level  of  232  dB  re  pPa/m  was 
used,  and  bathymetry  was  extracted  along  the  bearing  from  source  to  receiver  in  order  to 
accurately  simulate  the  true  environment.  The  sound  speed  profile  in  the  water  column 
used  by  the  model  was  taken  from  one  of  the  CTD  casts  in  the  region,  and  the  bottom 
was  assumed  to  be  homogenous  with  a  sound  speed,  density,  and  attenuation  of  1800  m/s, 
2.0  g/cm^,  and  0.2  dB/m/kHz  respectively. 

For  Segments  016, 031,  and  046,  the  near-bottom  vertical  line  array  (VLA)  was 
stationary  as  the  Cory’s  source  moved  around  it.  The  UTM  coordinates  (zone  23,  based 
on  the  WGS-84  ellipsoid)  for  the  near-bottom  array  were  (377.43,  2888.51).  The  source 
for  Segment  016  was  located  at  UTM  coordinates  (383.82,  2901.80)  so  that  the  range 
from  source  to  receiver  was  14.75  km.  The  source  for  Segment  031  was  located  at 
(398.03,  2899.89)  and  the  separation  range  was  23.53  km.  The  source  for  Segment  046 
was  located  at  (388.66,  2885.01)  and  the  separation  range  was  11.76  km. 

For  Segments  544  and  564,  the  near-bottom  VLA  was  located  at  UTM  coordinates 
(192.27,  2944.86).  The  source  for  Segment  544  was  located  at  (223.2,  2953.54)  and  the 
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separation  range  was  32.12  km.  The  source  for  Segment  564  was  located  at  (210.36, 
2932,48)  and  the  separation  range  was  21.92  km 

Plots  of  the  transmission  loss  fields  of  the  ocean  cross-sections  along  the  radials 
of  the  five  segments  were  made  at  the  center  frequency  of  227.5  Hz.  These  fields  are 
calculated  by  the  UMPE  model  and  show  the  propagation  of  the  acoustic  energy  from  the 
source  outward. 

Segment  016’s  arrangement,  shown  in  Fig.  3,  places  the  near-bottom  array  at  14.75 
km,  which  is  prior  to  one-half  convergence  zone  (CZ)  distance.  The  bottom  profile  for 
this  radial  is  relatively  smooth  with  a  slight  upslope  along  the  latter  half  of  the 
calculation  distance. 

Segment  03rs  arrangement,  shown  in  Fig.  4,  places  the  near-bottom  array  at  23.53 
km,  just  prior  to  one-half  CZ  distance.  The  propagation  in  this  case  is  bottom-limited  at 
25  km.  The  bottom  profile  for  this  radial  shows  an  almost  flat  bottom  along  the  mid¬ 
point  with  gradual  rises  of  approximately  1500  meters  near  source  and  receiver  locations. 
The  forward  scatter  resulting  from  interaction  with  the  bottom  can  be  dimly  seen  in  the 
top  right  of  the  figure  of  which  some  reflects  off  the  surface  back  down  to  the  location 
of  the  array  near  the  bottom. 

Segment  046’s  arrangement,  shown  in  Fig.  5,  places  the  near-bottom  array  at 
11.76  km  which  is  well  before  the  one-half  CZ  range.  The  bottom  profile  shows  several 
larger-scale  bottom  features  in  the  first  half  of  the  radial  distance. 

Segment  564’s  arrangement,  shown  in  Fig.  6,  places  the  receiver  array  at  21.92 
km,  prior  to  one-half  CZ  distance.  The  bottom  profile  shows  a  gradual  downslope  from 
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Figure  3.  Segment  016  Transmission  Loss  Field  for  the  Array  Source. 
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Figure  4.  Segment  031  Transmission  Loss  Field  for  the  Array  Somce. 
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Figure  5.  Segment  046  Transmission  Loss  Field  for  the  Array  Source. 
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Figure  6.  Segment  564  Transmission  Loss  Field  for  the  Array  Source. 


5-10  km  and  a  gradual  upslope  over  the  last  quarter  of  the  radial  distance,  with  no  large- 
scale  features  interrupting  the  relatively  smooth  profile. 

Segment  544’s  arrangement,  shown  in  Fig.  7,  has  one  of  the  most  interesting 
bottom  profiles  with  several  major  large-scale  features  from  which  reflected  energy  paths 
can  be  clearly  seen.  The  receiver  array  is  located  at  32.12  km,  just  past  one-half  CZ 
distance.  In  this  case  propagation  is  bottom-limited  at  -30.5  km,  prior  to  the  array 
location.  Figure  8  shows  the  UMPE  model  prediction  of  the  approximately  200-255  Hz 
broadband  pulse  arrival  transmission  loss  structure  over  the  entire  water  column  at  the 
range  of  the  near-bottom  VLA.  The  separation  between  primary  and  secondary  bottom 
interactions  is  observed  to  be  roughly  1.5  seconds. 

B.  EXISTENCE  OF  SECONDARY  MULTIPLE  ARRIVALS 

Imtially,  the  transmitted  pulse  in  all  cases  can  be  considered  a  single  wavefront 
in  time.  The  highest  energy  level  is  transmitted  horizontally  from  the  source  and  occupies 
the  higher,  more  shallow  end  of  the  pulse.  As  the  pulse  propagates  the  bottom  end 
undergoes  a  bottom  reflection  which  begins  to  spread  the  signal’s  range  resolution. 
Further  dowmrange,  this  forward  scattered  energy  reflects  from  the  surface  and  again 
interacts  with  the  bottom.  This  propagation  is  explained  in  detail  in  Ref.  [4]. 

The  range  resolution  of  this  secondary  interaction  has  been  much  degraded  as  the 
pulse  has  become  spread  out  over  some  range.  Also,  this  secondary  interaction  is 
occurring  at  the  same  time  the  primary  direct-path  pulse  is  ensonifying  a  patch  of  seafloor 
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Figure  7.  Segment  544  Transmission  Loss  Field  for  the  Array  Source. 
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Figure  8.  Segment  544  Broadband  Arrival  Structure  at  32.12  km. 


further  downrange.  Since  the  reciprocal  paths  of  both  the  primary  and  secondary 
interactions  back  to  a  monostatic  receiver  platform  will  arrive  at  the  same  time,  discrete 
scattering  events  at  the  position  of  a  primary  interaction  will  be  contaminated  or 
combined  with  scattering  events  at  different  places  on  the  seafloor.  These  secondary 
interactions  occur  at  higher  angles  of  incidence  and,  because  of  this,  may  have  their 
scattering  strengths  enhanced. 

The  effects  of  these  secondary  interactions  are  most  significant  at  ranges  of  one- 
half  CZ  and  greater.  In  cases  where  the  propagation  is  not  bottom-limited,  or  where  the 
primary  portion  of  the  pulse  is  forward  scattered,  the  highest  level  of  ensonification 
beyond  the  range  of  primary  interactions  will  be  due  to  the  secondary  interactions  instead 
of  the  primary.  The  returns  from  these  interactions  can  result  in  reverberation  features 
being  mapped  beyond  the  true  point  of  contact  by  several  kilometers. 

By  looking  at  the  data  recorded  by  the  near-bottom  VLA’s  in  the  vicinity  of  one- 
half  CZ,  such  phenomena  can  be  verified.  The  ability  to  beamfotm  the  data  allows  the 
primary,  near-grazing  interaction  to  be  distinguished  from  the  secondary,  higher-angle 
arrivals  undergoing  multiple  reflections. 

Data  from  the  five  simulated  segment  runs  has  been  extracted  over  the  aperture 
of  the  actual  near-bottom  VLA  (120  meters  long,  approximately  100  meters  off  of  the 
bottom).  Vertical  spacing  of  the  model  data  required  63  depth  elements  to  be  extracted 
to  simulate  the  complete  aperture  of  the  near-bottom  VLA.  This  data  was  then 
beamformed  as  explained  in  Chapter  II  to  yield  plots  of  transmission  loss  arrival 
structures  for  the  five  segments. 
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Figure  9  shows  the  beamformed  predicted  arrival  structure  for  Segment  016. 
Values  of  transmission  loss  in  dB  (unnormalized)  outside  the  plotted  scale  are  assigned 
the  maximiun  or  minimum  values  for  the  plot  as  appropriate.  The  exact  normalization 
factors  for  transmission  loss  (TL)  values  are  unknown  due  to  the  process  of  beamforming 
the  data,  but  all  modeled  and  beamformed  values  are  accurate  relative  to  each  other. 
Positive  angles  indicate  propagation  from  the  bottom  towards  the  surface,  i.e.  energy 
coming  up  from  below  these  element  depths.  The  initial  direct  path  arrival  is  so  intense 
that  the  sidelobe  structure  has  spread  over  multiple  adjacent  angles  and  times.  The  next 
set  of  arrivals  occurs  roughly  2.5  seconds  later  and  exhibits  a  symmetry  about  grazing 
with  lower  intensity  and  more  diffuse  arrivals  from  below  (positive  angles).  The  arrivals 
from  below  are  more  diffuse  due  to  an  extra  reflection  from  the  bottom  just  prior  to 
reaching  the  array.  These  secondary  arrivals  tend  toward  higher  arrival  angles  with 
increasing  time  as  expected  due  to  the  increased  path  lengths  they  must  travel.  In 
addition  to  these  features,  there  is  also  some  structure  which  tends  toward  more  grazing 
angles  with  time.  This  seems  to  originate  from  a  point  of  high-intensity  incidence  which 
also  produces  the  expected  increasing  angle  with  time  features.  The  cause  of  this 
unexpected  feature  is  not  well  understood,  though  all  features  must  be  the  result  of 
forward  scatter  since  the  model  does  not  compute  any  back-scattered  effects. 

The  match-filtered  and  time-domain  beamformed  LFM  signal  arrival  structure 
received  by  the  VLA  for  Segment  016  is  shown  in  Fig.  10.  The  primary  arrival  occurs 
at  0.4  seconds  and  shows  earliest  arrivals  from  above  (negative  angles)  with  later  arrivals 
occurring  as  the  arrival  angle  passes  through  zero  to  positive  values,  which  also  matches 
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Angle  From  Horizontal  Angle  From  Horizontal 


with  model  predictions.  The  structure  seen  at  10*  and  1.3  seconds  is  most  likely 
backscatter  from  some  bottom  feature  just  past  the  VLA.  This  structure  does  not  show 
up  in  Fig.  9  as  the  model  only  calculates  forward  scatter.  Secondary  arrivals  occur 
approximately  2.8  seconds  after  the  primary  and  exhibit  the  same  structure  as  described 
above  for  the  model  predictions  of  Fig.  9. 

Figure  1 1  shows  the  beamformed  predicted  arrival  structure  for  Segment  031.  The 
general  pattern  is  the  same,  though  the  second  set  of  arrivals,  which  occurs  roughly  1.5 
seconds  after  the  primary  arrival,  is  more  spread  out  in  time  as  would  be  expected  for  the 
longer  propagation  range. 

The  match-filtered  and  time-domain  beamformed  LFM  signal  arrival  structure 
received  by  the  VLA  for  Segment  031  is  shown  in  Fig.  12.  The  primary  arrival  occurs 
at  0.25  seconds,  there  exists  some  backscatter  structure  at  1.5  seconds,  and  secondary 
arrivals  show  up  1.8  seconds  after  the  primary  arrivals.  This  shows  good  overall 
agreement  with  the  model  predictions  shown  in  Fig.  11. 

Figure  13  shows  the  beamformed  predicted  arrival  structure  for  Segment  046. 
Secondary  arrivals  first  appear  2.8  seconds  after  the  primary.  The  match-filtered  and 
time-domain  beamformed  LFM  signal  arrival  structure  received  by  the  VLA  for  Segment 
046  is  shown  in  Fig.  14.  Again,  arrival  angles  and  time  between  primary  and  secondary 
arrivals  agree  very  well. 

Figures  15  and  16  show  the  beamformed  predicted  arrival  structme  for  Segment 
544.  No  actual  beamformed  data  was  available  for  this  segment.  The  model  predictions 
conducted  have  a  frequency  spacing  which  results  in  a  total  time  window  of  15.36 
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Figure  16.  Segment  544  Complete  Beamformed  Arrival  Structure. 
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seconds,  and  this  longer  time  window  is  shown  in  Fig.  16.  There  exists  a  third,  fourth, 
and  possibly  a  fifth  distinct  arrival  structure,  each  of  which  arrives  at  progressively  higher 
angles  and  with  less  time  separating  each  subsequent  set  of  arrivals.  Comparison  of  the 
full  time  windows  of  beamfotmed  modeled  data  shows  that  in  general,  the  farther  out  in 
range  from  the  source,  the  shorter  the  time  interval  between  the  primary  and  first 
secondary  arrivals;  the  more  spread  out  in  time  these  secondary  arrivals;  and  the  shorter 
the  time  between  subsequent  sets  of  secondary  arrivals,  all  of  which  must  result  from 
forward  scatter.  For  the  purpose  of  this  thesis,  only  the  secondary  multi-path  arrivals 
which  occur  at  or  close  to  the  same  time  as  the  primary  arrival  are  of  concern. 

Figure  17  shows  the  beamformed  predicted  arrival  structure  for  Segment  564  and 
Fig.  18  shows  the  match-filtered  and  beamformed  data  from  the  near-bottom  VLA  for  this 
segment.  Given  that  the  predicted  arrival  structures  for  all  five  segments  are  the  result 
of  only  forward  propagation  calculations,  there  is  significant  agreement  between  predicted 
and  actual  arrival  structures  for  all  cases,  showing  that  these  secondary  arrivals  do  indeed 
exist  and  can  affect  two-way  mapping  resolution. 

C.  INFLUENCE  OF  SECONDARY  INTERACTIONS 

The  existence  of  secondary  arrivals  has  been  clearly  shown  in  the  previous  section. 
To  what  extent  do  these  secondary  arrivals  degrade  the  return  signal  at  the  receiver  in  a 
monostatic  arrangement? 
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In  an  attempt  to  answer  this,  exact  two-way  travel  time  arrival  signal  levels  were 
calculated  for  Segments  016  and  544  as  explained  in  Chapter  11.  In  order  to  ensure  that 
a  Fourier-transformed  time  window  would  be  large  enough  to  capture  the  complete 
structure  of  the  two-way  propagation  out  to  the  farthest  range  needed,  frequency  spacing 
for  the  UMPE  modeling  runs  to  determine  exact  two-way  travel  times  was  halved  so  that 
there  were  2048  frequencies  centered  at  227.5  Hz  with  a  total  bandwidth  of  66.66  Hz, 
resulting  in  a  total  Fourier-transformed  time  window  of  30.72  seconds. 

The  exact  two-way  travel  time  structures  and  corresponding  bottom  profiles  for 
Segments  016  and  544  are  shown  in  Figs.  19-22.  There  is  much  more  structure  in  both 
Figs.  19  and  21  than  exists  for  the  corresponding  bottom  profiles.  TTus  is  an  indication 
of  the  high  degree  of  complexity  in  the  ensonifymg  field.  The  fact  that  such  complexity 
can  exist  distinct  from  any  fine-scale  bathymetric  roughness  is  described  in  more  detail 
in  Ref.  [4].  Large  features  can  be  discerned  by  large-scale  arrival  structures  and  trends 
in  the  return  signal,  but  small-scale  features  corresponding  to  the  "spilsy-ness"  of  the 
return  signal  simply  do  not  exist  in  the  modeled  bathymetry.  Lobe  structure  from  the 
source  array  as  seen  in  Fig.  3  is  the  cause  of  the  return  signal’s  general  appearance  of 
very  low  return  up  to  approximately  4  km,  stronger  returns  from  4  to  10  km,  and 
strongest  returns  from  10  to  15  km,  the  maximum  range  considered. 

In  an  attempt  to  remove  the  majority  of  secondary  interactions  from  the 
propagation,  the  UMPE  code  was  modified  for  each  segment  to  effectively  remove  the 
bottom  interface  up  to  a  range  of  five  kilometers  in  front  of  the  near-bottom  VLA’s.  This 
was  done  by  matching  water  and  bottom  densities  and  sound  speeds  along  the  radial 
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range.  At  five  kilometers  prior  to  the  VLA,  the  bottom  properties  were  linearly  changed 
so  that  within  three  kilometers  of  the  VLA,  the  bottom  properties  were  the  normal  values 
used  in  the  previous  runs  as  stated  earlier  in  the  chapter.  Though  the  bottom  interface 
was  removed,  energy  entering  the  ocean  bottom  was  still  attenuated  as  before.  The 
results  of  altering  the  bottom  properties  for  these  segments  on  the  propagation  of  acoustic 
energy  can  be  seen  in  Figs.  23  and  24  which  show  the  transmission  loss  fields  of  the 
ocean  cross-sections  along  the  radials  of  the  two  segments  at  the  center  frequency  of 
227.5  Hz. 

Figure  23  shows  the  propagation  along  Segment  016.  Surface  reflections  still 
occur,  but  any  energy  reflected  from  the  surface  near  the  source  would  hit  the  bottom 
interface  prior  to  the  VLA  in  the  altered  bottom  region  and  be  attenuated.  With  multipath 
energy  removed,  the  propagation  paths  of  the  primary  energy  close  to  the  surface  can  now 
be  clearly  seen.  The  three  primary  propagating  lobes  of  the  vertical  array  source  are 
distinct,  and  the  interference  pattern  resulting  from  Lloyd’s  mirror  effect  of  the  surface- 
reflected  energy  within  the  lobes  is  also  more  distinct.  Figure  24  shows  the  propagation 
along  Segment  544.  Bottom-reflected  energy  just  prior  to  the  VLA  can  be  seen  to  reflect 
from  the  surface  and  travel  back  down  to  ensonify  the  region  beyond  31  km  range.  This 
region  sees  little  of  the  primary  energy  which  is  bottom-limited  around  31  km. 

Figures  25  and  26  compare  the  unaltered  bottom  signal  return  and  the  altered 
bottom  signal  return  for  Segment  016.  At  ranges  approaching  15  km,  the  same  large 
structure  exists  in  the  returns,  while  in  Fig.  25  it  can  be  seen  that  multipath  effects  tend 
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Figure  23.  Segment  016  Transmission  Loss  Profile  -  Altered  Bottom. 
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Figure  24.  Segment  544  Transmission  Loss  Profile  -  Altered  Bottom. 
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to  fill  in  the  areas  of  low  primary  signal  return,  adding  up  to  5  dB  of  return  m  some 
cases.  While  the  effect  of  the  secondary,  multipath  arrivals  is  not  very  significant  at  these 
short  ranges,  it  is  important  to  note  the  considerable  differences  in  details  m  the  region 
of  the  altered  bottom.  The  distinction  can  be  considered  the  difference  between 
computing  the  total  (incident,  reflected,  transmitted,  and  forward  scattered)  field  at  the 
interface  versus  computing  only  the  incident  field.  To  correctly  predict  reverberation 

levels,  the  total  field  at  the  interface  must  be  computed. 

Figures  27  and  28  compare  the  unaltered  bottom  signal  return  and  the  altered 

bottom  signal  return  for  Segment  544.  Again,  the  general  effects  of  source  array  lobe 
structure  can  be  seen  when  compared  to  Fig.  7.  In  the  range  of  approximately  17  to  21 
km,  signal  returns  in  the  unaltered  bottom  case  are  as  much  as  25  dB  higher  than  the 
returns  caused  by  primary  interactions  alone.  By  referring  back  to  the  CW  field  in  Fig. 
7,  it  is  noted  that  a  ridge  structure  is  near  16  km.  This  would  produce  a  shadow  region 
in  the  17-20  km  range  if  only  direct  path  ensonification  were  considered.  However, 
secondary  multipaths  are  predicted  to  significantly  increase  the  levels  in  this  region. 
Another  interesting  feature  is  that  the  primary  return  in  the  unaltered  bottom  case  at  the 
range  of  approximately  21  to  23  km  is  actually  lessened  by  the  multipath  effects.  This 
is  presumably  due  to  the  difference  between  computing  the  total  field  versus  only  the 
incident  field.  Once  again,  the  lower  primary  return  at  the  range  of  25  km  due  to  a 
shadow  zone  behind  a  ridge  at  approximately  23  km  has  had  as  much  as  30  dB  added  to 
it  in  Fig.  27  from  the  returns  of  secondary,  multipath  interactions  being  received  at  the 

same  time  as  the  primary  return. 
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Range  in  Kilometers 

Figure  28.  Segment  544  Two-Way  Exact  Travel  Time  Structure  -  Altered  Bottom. 


Figures  29  and  30  show  a  magnified  view  of  the  imaltered  bottom  signal  return 
and  the  altered  bottom  signal  return,  respectively,  for  Segment  544  for  the  ranges  of  25 
to  33  km.  The  secondary  signal  returns  in  Fig.  29  could  result  in  some  features  being 
mapped  to  incorrect  seafloor  locations.  The  structure  in  Fig.  29  from  31  km  to  32  km 
which  exhibits  the  gradual  upslope  then  downslope  is  actually  from  secondary  returns 
occurring  approximately  1.5  seconds  later  (see  Figs.  8  and  15)  from  the  large-scale  ridge 
face  at  30  km  and  would  result  in  mapping  another  ridge  structure  at  this  location.  This 
secondary  "echo"  was  previously  predicted  in  Ref.  [4]. 
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IV.  SUMMARY 


It  has  been  shown  that  the  secondary  arrivals  predicted  by  the  UMPE  model  do 
exist  in  the  measured  data.  The  exact  two-way  travel  time  calculations  show  that  the 
levels  of  ensonification  resulting  from  these  secondary  multipaths  can  be  significant 
compared  to  the  primary  direct-path  interaction  at  ranges  near  one-half  GZ  and  greater. 
Furthermore,  the  exact  two-way  levels  display  a  highly  complex  nature  which  are  due  to 
complicated  ensonification  structures  rather  than  small-scale  bathymetric  features.  The 
overall  implication  of  this  is  that  there  are  limitations  on  the  ability  to  resolve  discrete 
structures  at  ranges  beyond  a  few  water  depths.  These  limitations  arise  from  system 
limitations,  processing  techniques,  and  the  inherently  complex  nature  of  the  propagation. 

One  way  to  overcome  some  of  these  limitations  in  mapping  resolution  might  be 
to  use  a  vertical  receive  array,  or  some  form  of  volumetric  array,  so  that  the  ability  to 
beamform  in  the  vertical  exists.  Most  secondary  multipath  interactions,  at  least  those 
occurring  during  the  return  of  the  acoustic  signal,  will  be  arriving  at  higher  angles  than 
the  direct-path  returns  so  that  they  can  be  selectively  rejected  from  the  signal  return. 

Another  aspect  of  the  propagation  which  may  affect  the  resolution  of  such  a 
system  is  coupling  between  radials.  This  has  not  been  considered  in  this  research. 
Further  studies  incorporating  azimuthal  coupling  should  be  conducted  in  an  attempt  to 
better  understand  the  inherent  problems  associated  with  attempting  to  precisely  map 
bottom  bathymetry  in  this  manner. 


47 


Finally,  the  shallow  water  environment  is  much  different  from  the  deep  ocean 
environment  and  the  modeling  done  for  this  thesis  could  be  applied  to  shallow  water 
cases  to  determine  the  effects  of  these  secondary  interactions.  In  a  shallow  environment, 
most  of  the  propagation  will  be  bottom-limited  so  multipaths  will  be  the  primary  source 
of  bottom  ensonification.  With  the  shift  in  emphasis  toward  the  littoral  environment  for 
the  Navy  of  the  future,  these  issues  will  need  to  be  addressed  for  improved  sonar  system 
performance  on  the  submarines  and  ships  of  tomorrow’s  Navy. 
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